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Abstract 

The cumulants of thermal variables are of general interest in physics due to their extensivity and 
their correspondence with susceptibilities. They become especially significant near critical points 
of phase transitions where they diverge along with the correlation length. Cumulant measurements 
have been used extensively within the held of heavy-ion physics, principally as tools in the search 
for a hypothetical QCD critical point along the transition between hadronic matter and QGP. 
The volume of individual heavy-ion collisions can be only partially constrained and, as a result, 
cumulant measurements are significantly biased by the limited volume resolution. We propose 
a class of moments called strongly intensive cumulants which can be accurately measured in the 
presence of unconstrained volume fluctuations. Additionally, they share the same direct relationship 
with susceptibilities as cumulants in many cases. 
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I. INTRODUCTION 


The development of classes of moments that are either invariant or have clear scaling 
properties under various operations has been an active area of research for well over a 
century. The initial work on cumulants themselves predated the development of statistical 
partition functions and focused on their scaling properties under affine transformations [jl]. 
They were originally called half-invariants due to these scaling properties and were later 
renamed cumulants due to their additivity under convolution [2|. In the field of computer 
vision, two-dimensional moments which are invariant under operations such as rotation, 
scaling, and translation have played a central role in pattern recognition M- Within the 
field of particle physics, the R fluctuation measures were constructed as ratios of factorial 
moments such that detector efficiencies would cancel resulting in moments invariant under 
binomial efficiency losses 13 El- 

More recently, the A and £ observables were proposed to address the issue of measurement 
biases resulting from the poor constraint of volume in heavy-ion collisions |ZHS] • These 
quantities are two-dimensional second-order moments which are constructed in such a way 
that the volume fluctuation terms cancel. The authors coined the term strongly intensive 
to describe these observables because they are not only independent of the volume of the 
system but also of the distribution of volume within an ensemble. They have been used 
effectively in fluctuation analyses [lU], 11], but their physical meanings are obscure relative 
to those of cumulants and they only measure second-order fluctuations. 

In this work, we present a new set of statistical quantities which we call strongly in¬ 
tensive cumulants. These quantities are invariant under both convolution and mixing with 
distributions sharing the same strongly intensive cumulants. A direct result of this is that 
they can accurately be measured experimentally in situations where the volume is not well 
constrained. These quantities are directly related to cumulants when certain conditions are 
met and therefore, in physical systems, to thermodynamic susceptibilities. 

The strongly intensive cumulants are of particular interest in heavy-ion collisions where 
cumulant measurements of conserved charges have been proposed as a signature for critical 
point fluctuations m The experimental determination of these cumulants has been a focus 
of the RHIC beam energy scan dang and the published measurements have been discussed 
in the context of critical fluctuations as well as compared to lattice calculations in order to 


2 


determine the temperature and chemical potentials at chemical freeze-out p5|. It has been 
shown that these measurements depend greatly on the method used to constrain collision 
volumes which calls into question the validity of any physics conclusions drawn from the 
results pH]. 

To demonstrate the utility and efficacy of the strongly intensive cumulants we present 
an analysis of net proton fluctuations in simulated heavy-ion collisions. This analysis makes 
clear the issues with current analysis techniques and illustrates how the strongly intensive 
cumulants resolve them. In light of this, we propose that the strongly intensive cumulants 
listed in Eq. ([9]) be used as drop-in replacements in future heavy-ion cumulant analyses. 

II. MOTIVATION 

We begin by briefly reviewing the generating function formalism with respect to sta¬ 
tistical moments, cumulants, and their relationships to volume fluctuations. Let X = 
{X \, X 2] ... ,X n ) be the components of a random vector. The moment-generating function 
is then defined to be 


<M0 = ( e&Xi ) x 
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where the (• • ■) x notation denotes the expectation value over the distribution of X and the 
repeated index i is implicitly summed over according to Einstein summation notation. The 
coefficients of the Taylor series correspond to the moments of the distribution and can be 
recovered from the generating function by taking derivatives and setting £ = 0. Defining 
the operator D, to be a /d£ il we can see that 


Hn.r, ., r „ = K'D? Wle,0 

= (Xl'X?-X?) x 


as we would expect. This approach of recovering coefficients from the generating function 
can be applied in the same way for both the standard cumulants and the strongly intensive 
cumulants that we will introduce later. 
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The generating function for the cumulants is then defined in terms of the moment¬ 
generating function as 


0(f) = In 0(0 
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where the coefficients Ki,j,...,k are the cumulants and can be expressed in terms of the moments 
by matching terms in the Taylor series, either combinatorially or using recursion relations 

HU- 


The utility of the cumulants becomes clear when we consider the convolution of probabil¬ 
ity distributions. If we take two independent random vectors, X and Y, then the distribution 
of Z = X + Y is given by the convolution of their respective distributions. The moment¬ 
generating function of Z is then given by <fiz (£) — { e ^ X+Y ^)x y w hich, due to the indepen¬ 
dence of X and Y, can be factored as (f>z (£) = ( e ^ iXi ) x (e^ Yi ) y . The cumulant-generating 
function is then simply 'ijiz (f) = In (pz (£) = In <j> x (£) + In 0y (£) which is precisely the sum 
of the cumulant-generating functions for A" and Y. Thus, we see that the cumulants of 
two probability distributions are additive under convolution, a fact that led to their current 
name. 

This additivity property of cumulants is closely related to their utility in physics. If 
we consider two volumes of matter that are each in the same thermodynamic state then 
the distributions of any total quantities (e.g. net charges, total energy) for the combined 
volume will be given by the convolution of their distributions for each of the two independent 
volumes. The cumulants of these distributions will necessarily be extensive and, after scaling 
by the volume, will give intrinsic quantities determined by the thermodynamic state of the 
matter. By relating the partition function to a moment-generating function we can see 
that the coefficients of the corresponding cumulant-generating function then encode how 
the mean values of total quantities change with respect to state variables like energy density 
or chemical potential [HB] . 

The determination of these coefficients is of general interest in physics but there is an 
experimental limitation that makes it difficult. If the volume of an ensemble of systems 
cannot be perfectly constrained then the cumulants of the distributions depend strongly 
on the distribution over volume within the ensemble. To see this we define an intrinsic 
generating function such that 
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where the primes indicate independence from the volume. We use the term volume here 
loosely; it really corresponds to any measure with which the cumulants scale linearly. In 
actuality, V can depend on the temperature, the energy density, and other quantities in 
addition to the volume. The important point in the context of this discussion is simply that 
this quantity can be factored out from the cumulant-generating function, as done in Eq. ([TJ). 

The moment-generating function can then be expressed as 


<H0 = e V *' t(} 

= 1 + vy «) + iry K) 2 + ■ • • 

with an explicit volume dependence. In any experimental context, the moments that are 
measured are those of a mixed distribution over an ensemble of volumes. The measured 
cumulants will then be described by the generating function 

^(0 = ln(e' / *' K) ) v , (2) 

= In (l + { V) v V>' «) + i {V 2 ) v iP' (?) 2 + • • •) 

which can easily be seen to not equal {V) v i/)' (£) unless V is fixed at a single value. 

The measured cumulants instead depend on the distribution of the volume in a straight¬ 
forward way [ID]. This relationship can be made more clear by considering the simple 
example of the variance of a single variable. By expanding the logarithm in Eq. ([2]) and 
matching terms we End that the measured variance would be 


K2=(V) v k.' 2 +({V 2 ) v -(V)1) ( k ',) 2 


( 3 ) 


in the presence of volume fluctuations. If the volume is constrained to a single value then 
the variance of V goes to zero and the relationship reduces to k 2 = (V)y ^4, the quantity 
that one would truly want to measure. The measured variance will be artificially high for 
any other distribution of volume due to the contribution of the second term. The most 
extreme example of this is a situation where k' 2 = 0, which would be approximately true for 
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the total number of atoms in a crystal lattice. In this scenario, the first term on the right 
hand side of Eq. (|3]) vanishes and the measured variance would be directly proportional to 
the variance of the volume in the measurement ensemble. 

It may seem as though the fluctuation terms could be subtracted off from measurements 
such that the cumulants (V) v k[ could be directly determined. This is true in theory, but 
the measurements cannot be corrected for without an exact knowledge of the volume dis¬ 
tribution. The precise shape of this distribution is typically not known in practice and so 
this approach is not applicable. Instead, measurements tend to be made without any at¬ 
tempt at corrections and, as a result, they are biased in poorly understood ways due to the 
contributions from volume fluctuations. 


III. DERIVATION 


We introduce here a set of statistical quantities called the strongly intensive cumulants. 
The goal is to construct a set of non-trivial statistical quantities that can be measured 
in physical systems without any dependence on the volume distribution. Formally, this is 
equivalent to saying that the quantities are invariant under both convolution and mixing of 
distributions sharing the same strongly intensive cumulants. We will begin by defining the 
strongly intensive cumulants and then prove that they satisfy the desired properties. 

Their generating function, 0*, is defined in terms of the partial differential equation 


d u (r) 


Du ( 0 ) 
D v ( 0 ) 


( 4 ) 


where the choice of u and v determines the exact flavor of the generating function. We assume 
from here foreword that the u and v indices are 1 and n, respectively. The components of 
X can always be rearranged such that this is the case, with the trivial exception of u = v. 
Without loss of generality, we can assume this has been done. 

The choice of the first and last components of X is arbitrary, but we will see in Sec. [V] 
that it makes no difference which component comes first with an appropriate choice of X n . 
In this situation, X n serves as a measure of the volume and is independent of the quantities 
that one is primarily interested in. The general case is more subtle but thinking of a physical 
situation where X n is a noisy volume measurement and the choice of X\ is arbitrary can be 
helpful in understanding how we proceed. 
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The strongly intensive cumulants, fiVi,r 2 ,...,r 2 j then correspond to the Taylor series coeffi¬ 
cients of 



in the same way that cumulants are defined by the coefficients of ip. The Kq r2 r3 rn terms in 
this series are, as of yet, undetermined because they correspond to the integration constants 


obtained when integrating Eq. Q with respect to £ 1 . We choose to define these terms to be 


Kq Tn = k* 2 r 3 r n , the strongly intensive cumulant obtained when the first component 
of £ is dropped such that £ —» (£ 2 , £ 3 , • • • , £ n ). This process can be applied repeatedly until 
the first index is nonzero or, if all indices are zero, Kq 0 0 = 0. This definition allows us to 
prove properties for D\ (ip*) and have them apply more generally to ip* because the missing 
K o r 2 ,r 3 ,...,r„ terms correspond to coefficients of D\ (ip*) for a different choice of £. 

First, we’ll show explicitly that the strongly intensive cumulants are completely indepen¬ 
dent of the distribution over volume in a physical system. This can be proved very simply 
by substituting the volume mixed moment-generating function from earlier into Eq. (Jd]) 



and canceling the volume dependent terms. This proof of their strongly intensive property 


is straightforward, but it is not the most clear way to demonstrate why this happens. By 
exploring the properties of the strongly intensive cumulants under convolution and mixing 
it will become more clear. 

We will first consider convolutions, whereby a new random vector Z is constructed as 


the sum of two other random vectors X and Y. For physical systems, this operation can be 
viewed as constructing a larger volume out of two smaller volumes. Given that Z — X + Y, 


we find that ip z (£) = ipx (£) (0 or, equivalently, that (p z (£) = (px (£) 4>y (£)• Plugging 


this into Eq. Q we find that 



D n ( ip x ) D x (ip* x ) + D n (' ip Y ) D 1 (V>y) 


( 6 ) 
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which can be seen as a ^-dependent weighted average of the differential equations defining 
the strongly intensive cumulant generating functions of X and Y. 

Now let us consider the case that X and Y have the same set of strongly intensive 
cumulants. This implies that D\ (ip* x ) = D\ ( ip Y ) which simplifies Eq. ([6]) to 


_ Dn jlpx) D\ (ipx) + D n ('P’y) Dl ( 1 p * X ) 
D n (■ ipx) + D n (■ ipY) 

_ Dn (j>x) + Dn (ipy) n ( I* \ 

D n (iPx) + D n (ip Y ) Dl [ix) 

= Di (ip* x ) = Di (r Y ) 


(7) 


showing that the strongly intensive cumulants are invariant under the convolution of distribu¬ 
tions with identical strongly intensive cumulants. This demonstrates the intensive property 
of the strongly intensive cumulants and how it emerges from the way they combine under 
convolution. 

The situation for mixing distributions very closely parallels that for convolving them. If 
convolution can be thought of as an operation for constructing new volumes then distribution 
mixing can be thought of as the operation of combining different volumes into an ensemble. 
If we imagine that half of the time we choose a random vector Z according to Px (Z) and 
half of the time we choose it from Py (Z) then the probability for the mixed distribution 
Z is given by P z (Z) = \ (P x (Z) + P Y (Z)). We could have chosen a ratio other than 
7 but the extension to a mixing parameter is trivial. The resulting moment-generating 
function is, like the probability distribution, the arithmetic average of those for the separate 
distributions: (pz (0 = ^ (<Px (0 + <Py (0) • This results in a strongly intensive cumulant 
generating function of 


Di (t z ) 


D n ( 4 >x ) D i ( ipx ) + D n (( p Y ) Di (ip Y ) 
D n (4>x) + D n (<p Y ) 


( 8 ) 


which is again a weighted average of the two independent differential equations. 

It is important to note that Eq. (J8]) is identical to Eq. © with ip replaced by (p. Although 
the strongly intensive cumulants do not combine as simply as the moments under mixing or 
the cumulants under convolution, they have the same relationship with the moments under 
mixing as they do with the cumulants under convolution. Similarly, we find that 


Dl ((P) _ y^Di ((P) _ Di (In (ft) _ Dl (iP) 
D n (<p ) 1 /<t>D n (4>) D n (In (p) D n (ip ) 
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showing that the strongly intensive cnmnlants are directly related to the cumnlants and 
moments of a distribution in the same way (i.e. the expressions for the strongly intensive 
cumulants are unchanged when all moments are replaced with cumulants or vice versa). 
The similarity in the relationships that the strongly intensive cumulants have with both the 
moments and the cumulants is at the core of why they exhibit strongly intensive behavior. 

It follows from Eq.s ([6]) and ([8]) that the strongly intensive cumulants are invariant under 
mixing and convolution of distributions with identical strongly intensive cumulants. This 
was shown explicitly for convolution but the same proof employed in Eq. (J7]) applies to 
mixing when Eq. ([8]) is used as a starting point. The invariance under convolution implies 
that they will be intensive, or independent of volume, in a thermodynamic system. The 
combination of an ensemble with many different volumes is, in effect, distribution mixing 
and the strongly intensive cumulants will therefore be invariant under this operation as well. 
This is equivalent to saying that the strongly intensive cumulants can be measured over an 
ensemble of different volumes without any dependence on the volume distribution. This line 
of reasoning helps to illustrate why the volume terms were shown to cancel in Eq. ([5]). 


IV. OBTAINING EXPRESSIONS FOR THE STRONGLY INTENSIVE CUMU¬ 
LANTS 


We move now to the task of determining polynomial expressions for the strongly intensive 
cumulants. Let / (£) = D n ( 0 ) and g (£) = 1 // (£) where the coefficients in the Taylor series 
for g (£) are given by a ri! r 2 ,...,r n / (o!r 2 ! • • • rj). We then End that 
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when at least one of ri,r 2 ,..., r n are nonzero. We can rearrange this to give a recursion 
equation 


a ri ,r 2 , 


ri r 2 r n 
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with the starting point a 0 ,o,...,o = g (0) = 1 /mo,o,...,o,i- 

Now we can express the strongly intensive cumulants in terms of these coefficients in a 
similar manner. We find that 
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which provides a full prescription for finding the strongly intensive cumulants for rq > 0. If 
r'i = 0 then the hrst component can be removed, as described earlier. For the special case 
of n — 2 and rq = 0 this can be simplified to 


r— 1 


K r ~ K r ,0 — 


h 0.1 


hr,o y ^ 


i= 1 


h?—i,l 


where the second zero index has been dropped. This can be done without ambiguity because 
the strongly intensive cumulants are trivial in one dimension. This recursive equation can 
be used to fold the explicit expressions 


* hi,o 

«1 = - 
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* h 2,0 hl, 0 hl,l 

«2 = —-2“ 
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=--- 5 -- 1 - o - 
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* h4,0 3/^3 ohl,l + hl,0h3,l 3/12,0h2,1 

^4 o 9 

ho,i ho,i ho,i 

^_ 6 ^t 2 ,ohl,l 2 + 6 hl, 0 h 2 ,lhl,l _ 6 /lgohl,l 3 
ho,l 3 ho.I 1 

which we expect will be most practical in common usage. 


V. RELATION TO OTHER QUANTITIES 

The name “strongly intensive cumulants” may seem contradictory, as they are not cumu¬ 
lative like the standard cumulants. This name was chosen due to their close relationship 
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with cumulants. We’ll work with a fixed volume system so that there are no volume fluc¬ 
tuation terms anywhere. Now consider the case where X n is completely uncorrelated from 
X U X 2 , ..., x n _ x at hxed volume such that (X^X? ■ ■ ■ X r n n ~?) (X?) = (X[ l X r 2 2 ■ ■ ■ X£) 
for any choice of ri, r 2 , ■ ■ ■, r n . It then follows that 

Di ( 0 ) Di ( 0 ) Di ( 0 ) 


D i m = 


D n (<t>) D n ((eti*) x ) (X n e^) x 

Di ( 0 ) Di ( 0 ) 

{X n )x( e ^ Xl )x (^n).Y 0 
D\ (In 0) L>!(0) 


(Xn) 
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(Xn) 
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which directly implies 


k, 


K, 


ri,r 2 ,...r n 


r 1 ,r 2 ,...r n 

K 0 , 0 ,..., 0,1 


( 10 ) 


for r\ > 0 and n > 1. The definition of the integration constant was, however, chosen 
such that this equality holds for r\ = 0 as well. Thus, we find that the strongly intensive 
cumulants are equal to their corresponding cumulants normalized by a volume term given 
that X n is independent from the other components. For the case of cumulant ratios, which 
are used frequently in experimental contexts, this normalization cancels. 


When we extend this to a variable volume, we find that the left hand side of Eq. (10) 


does not change because it is strongly intensive. This means that it corresponds to what 
the cumulant ratio on the right hand side would be in the absence of volume fluctuations. 
In a very real sense, this quantity corresponds to what one would want to measure in the 
absence of volume fluctuations. 

It is important to note here that this correspondence with the cumulants depends only on 
the independence of X n at hxed volume and not at all on the shape of its distribution. This 
means that any estimate of the volume can be used for X n , regardless of its noise profile. 
Any independent quantity that has some dependence on the volume can be used as X n if it 
is adjusted such that the dependence is linearly proportional. This opens up the possibility 
of using volume estimates for that have been previously considered too noisy or irregular 
for volume determination. 

It is also worth mentioning that k* r2 r should equal zero when X n is chosen to be 
independent at hxed volume, given that at least one of ri,r 2 , ■ ■ ■, r n _i are nonzero. This 


follows trivially from the cumulant correspondence demonstrated in Eq. (10) because all joint 


cumulants involving X n must be zero in order for it to be entirely independent. This means 


11 












that in experimental contexts the strongly intensive cumulants with r n = 0 will likely be the 
most interesting. It also suggests that evaluating how far the strongly intensive cumulants 
with r n > 0 are from zero might be useful for evaluating systematic errors involving the 
independence of X n . 

The second-order strongly intensive cumulants are also closely related to the £ and A 
fluctuation measures that were mentioned earlier [MJ. Namely, 

£ PA, X 2 ] = k* 2 [X l5 X 2 ] + 4 [X 2 ,Xi] 

A [A 1; X 2 ] = k * 2 [X l5 Ad] - k * 2 [Ad, Ad] 

when the normalization factor on £ and A is chosen to be , and the [Xi,X 2 ] denotes 

the order of the random vector components. The strongly intensive cumulants have the 
same properties as these quantities but make the relationship to the underlying physics 
more explicit. 


VI. A PRACTICAL EXAMPLE 

We consider here a practical example from heavy-ion physics. Higher order fluctuations 
of conserved quantities in nuclear collisions has been a topic of considerable interest in recent 
years. The most discussed of these quantities has been the ratio of n 4/^2 for net proton 
fluctuations as measured by the STAR experiment. With respect to these measurements, 
it has frequently been reported that “there are interesting trends, including e.g. the drop 
in the kurtosis of the net-proton distribution at y/sNN= 27 and 19.6 GeV” (20]. We will 
study simulated nuclear collisions to quantify the magnitude of bias introduced by volume 
fluctuations and to demonstrate the utility of the strongly intensive cumulants for addressing 
this issue. 

Our analysis was run over 35 million Au+Au events at y/sjsrN = 7.7 GeV that were 
generated using the UrQMD model Ell 122]. The procedure followed was designed to closely 
mirror those used by STAR mm- The net proton number in each event was quantified 
as the number of protons minus the number of antiprotons with transverse momenta in 
the range 0.4 GeV < px < 0.8 GeV and pseudorapidities between —0.5 < r] < 0.5. The 
detector efficiency for these particles was assumed to be unity in order to disentangle volume 
fluctuations from other effects. The approximate volume of the system was measured using 
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the multiplicity of charges particles with 0.5 < [)]\ < 1.0 and px > 0.15 GeV, a quantity 
called Ref mu i t:2 within STAR. This is the centrality quantity used in the net charge, rather 
than net proton, analysis at STAR but it was chosen to allow for a measurable X n variable. 

Centrality was defined as the integrated percentile of Ref mu i t: 2 , from the largest values 
to the smallest. The cumulant ratios were computed for each individual value of Ref mu i t: 2 
and then averaged across each centrality bin. This was done to minimize the impact of 
binning on the results. This procedure was repeated three times with different binomial 
efficiencies for Ref mu i t: 2 each time: p = 1/3, p = 2 /3, and p — 1. The efficiency of p = 2 /z 
is the most realistic but the other two were included to demonstrate the effect of better 
and worse volume resolutions. Additionally, the analysis was repeated with a multiplicity 
variable that counted all pions and kaons produced in each event. The results obtained using 
this multiplicity variable represent an ideal case of resolution where the biases induced by 
volume fluctuations should be largely eliminated. 

Finally, the analysis was repeated using strongly intensive cumulants rather than standard 
cumulants. The X t component corresponded again to the net proton number while X 2 was 
chosen to be the number of charged pions and kaons with 0.4 GeV < px < 0.8 GeV and 
—0.5 < 7] < 0.5 . The results can be seen in Fig. [TJ 

The most striking feature of the results is how dramatically the cumulant measurements 
are shifted by volume fluctuations, in some cases by well over a factor of two. The features 
observed in the STAR data are on order of 20-40% which can easily fall within the unquan¬ 
tified systematic errors caused by volume fluctuations. This is particularly true considering 
that t/snn = 27 and 19.6 GeV were run with different detector configurations than the other 
Beam Energy Scan (BES) energies and therefore have significantly different systematics. 

Another interesting thing to note is the drastic shift in the cumulant ratios in the most 
central events. This same trend can be seen in the data and it can be explained by a 
suppression of volume fluctuations when selecting on the highest multiplicity events. For 
mid-central events, a single multiplicity value corresponds to high multiplicity fluctuations 
from smaller volumes and low multiplicity fluctuations from larger volumes. This results in 
a wide range of possible volumes. For the highest multiplicity values there are only upward 
fluctuations due to there being roughly a maximum volume attainable in a collision. This 
leads to a tighter constraint of the volume near this maximum value for the most central 
events. This, at first glance, might suggest that the most central collisions are the most 


13 




Figure 1. Net Proton Cumulant Ratios 

Standard and strongly intensive cumulant ratios of the net proton distribution in UrQMD 
generated Au+Au events at ^/snn = 7.7 GeV as a function of centrality. The square markers 
indicate standard cumulant ratios and the circle markers indicate the strongly intensive cumulant 
ratios for various centrality resolutions. The points are offset slightly along the x —axis for clarity. 
The blue band indicates what the standard cumulant ratios would be with an ideal centrality 
resolution. We see that the strongly intensive cumulants correspond very closely with the ideal 
centrality resolution scenario regardless of the actual centrality resolution. The standard cumulant 
ratios, on the other hand, are very significantly biased and depend quite strongly on the centrality 
resolution. This shows very clearly that the strongly intensive cumulants more accurately measure 
the desired cumulant ratios in the inevitable presence of volume fluctuations. 

trustworthy. Unfortunately, this is not the case in practice because these same bins are the 
most sensitive to pileup events and secondary collision background. 

Moving on to the strongly intensive cumulants we find that they are nearly identical for 
each of the resolution settings. Some small differences can be seen due to physics correlations 
between the protons in X\ and the pions and kaons in X 2 . but, even still, the effects of volume 
fluctuations are suppressed by well over an order of magnitude compared to the standard 
cumulant ratios. The remaining differences could be largely eliminated by a more careful 
choice of X 2 . These quantities aren’t only invariant, they can also be seen to match almost 
exactly with what the cumulant ratios would be in the case of ideal centrality resolution. 
This very clearly illustrates both the invariance under volume fluctuations of the strongly 
intensive cumulants as well as their correspondence to the cumulants in the absence of 
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volume fluctuations. 


VII. CONCLUSIONS 

A new set of statistical moments called the strongly intensive cumulants have been pre¬ 
sented. These quantities have been shown to be invariant under both convolution and mixing 
with distributions sharing identical strongly intensive cumulants. This property allows them 
to be experimentally determined in ensembles of physical systems where the distribution over 
volume is unknown. We have studied a practical example from heavy-ion physics where the 
measurement of cumulant ratios has been shown to be extremely biased due to the uncer¬ 
tainties in constraining the system volume. In this same example, the strongly intensive 
cumulant ratios have been shown to be almost entirely independent of how well constrained 
the volume is. Furthermore, they have been shown to correspond to the cumulant values 
that would be measured if the volume could be almost perfectly constrained. For these 
reasons, it is clear that the strongly intensive cumulants are better suited for determining 
underlying physics in systems with imperfectly constrained volumes. 
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